Exchange coupling in Eu monochalcogenides from first principles 
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Using a density functional method with exphcit account for strong Coulomb repulsion within the 
4/ shell, we calculate effective exchange parameters and the corresponding ordering temperatures of 
the (ferro) magnetic insulating Eu monochalcogenides (EuX; X=0,S,Se,Te) at ambient and elevated 
pressure conditions. Our results provide quantitative account of the many-fold increase of the 
Curie temperatures with applied pressure and reproduce well the enhancement of the tendency 
toward ferromagnetic ordering across the series from telluride to oxide, including the crossover from 
antiferromagnetic to ferromagnetic ordering under pressure in EuTe and EuSe. The first and second 
neighbor effective exchange are shown to follow different functional dependencies. Finally, model 
calculations indicate a significant contribution of virtual processes involving the unoccupied / states 
to the effective exchange. 

PACS numbers: 75.50.Pp,75.30.Et,71.27.+a 



Ferromagnetic semiconductors have become object of a 
great technological interest with the appearance of spin- 
tronics because they can provide a spin-dependent tun- 
neling barrier. Especially challenging is to achieve a siz- 
able ordered moment at room temperature, which is cru- 
cial for a large scale application of the technology. The 
ability to calculate the ordering temperature and under- 
stand the exchange mechanisms on the material specific 
level is of particular importance. There is currently an 
intense effort to locate such materials within the dilute 
magnetic semiconductors, where the magnetic moment 
is carried by impurities in an otherwise non-magnetic 
system, but few candidates have been found. An at- 
tractive alternative is EM insulators. Ferromagnetism is 
rare in stoichiometric materials without charge carriers. 
Europium monochalcogenides (EuO, EuS, EuSe, EuTe) 
belong to this small group of ferromagnetic insulators [l| . 

Crystallizing in the rock-salt structure, the first two 
members of the group order ferromagnetically at 69.2 K 
and 16.6 K respectively [3, while EuTe, a type II antifer- 
romamet, becomes ferromagnetic only at elevated pres- 
sure Q . EuSe is at the border line between ferromagnetic 
and antiferromagnetic order with ferromagnetism stabi- 
lized by a moderate pressure of 0.5 GPa ^. Application 
of pressure strongly enhances Curie temperatures of all 
these materials. The Eu^+ valency results in half filling of 
the Eu 4/ shell with ^S configuration of the groundstate 
multiplet. Due to the localized nature of the moment- 
carrying / orbit als, inter-site exchange interactions can 
only be mediated by the valence and conduction elec- 
trons. The intra-atomic f — d and f — s exchange, which 
is the leading /-valence interaction, gives rise to temper- 
ature dependent features (red shift effect) in the valence 
electron spectrum which are well captured by the ferro- 
magnetic Kondo lattice model [5|, p. It is due to the 
insulating groundstate that this interaction alone cannot 



give rise to an effective coupling between the local mo- 
ments. To do so excitations across the gap and/or mixing 
of the / and valence/conduction bands has to be taken 
into account. The relevant exchange processes have been 
discussed previously on a qualitative level 7, 8, 9]. While 
a significant amount of ah initio calculations of the Curie 
temperature in metallic systems has been done (e.g. Ref. 
10 and references therein) attempts to address the Curie 
temperature and coupling mechanisms of ferromagnetic 
insulators on a first principles level are rare and become 
quite involved IL _12j. 

The / states in rare earths such as Eu pose an extra 
challenge in obtaining a quantitative microscopic theory. 
Ah initio electronic structure methods based on density 
functional theor y (D FT) _13 ] and the standard semi-local 
approximations [ij, |l5| have notorious problems in deal- 
ing with the strong correlations within the 4/ shell. In 
particular these approximations often result in incorrect 
filling of the 4/ states. Early bandstructure calculation of 
europium chalcogenides by Cho 16] using empirical po- 
tential did not address the inter-site exchange coupling. 
In this work we use the LDA+U method 17], which pro- 
vides correct filling of the 4/ states while allowing for 
mixing with the rest of the band (unlike the frequently 
used open-core treatment). Just as importantly, posi- 
tions of the resulting occupied and unoccupied / bands 
are realistic. The half filling of the 4/ shell in Eu^+ 
removes additional problems associated with orbital de- 
grees of freedom. 

The calculations reported here were performed using 
the Wien2k 18] implementation of the full-potential lin- 
earized augmented-plane-waves (FLAPW) method with 
the rotationally invariant LDA+U functional and double- 
counting scheme of Ref. [rj The size of APW+lo basis 
was determined by the cut-off RmtKmax=8 correspond- 
ing to approximately 100 basis functions per atom. Ap- 
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FIG. 1: The nearest-neighbor Ji and next- nearest-neighbor 
J2 effective exchange parameters in EuO and EuS as func- 
tions of the lattice constant calculated for different values of 
U (circle - 6 eV, square - 7 eV, diamond - 8 eV, triangle - 9 
eV, the lines serve as guides for eye). The deviation for the 
linear dependence of Ji at high pressures for EuO coincides 
with the onset of metallic behavior due to overlap of the / 
and valence band. 



proximately 30 irreducible k-points (depending on the 
magnetic structure) out of the 250-k-point regular grid 
were used in the Brillouin zone integrations. The cal- 
culations were performed for lattice constants spanning 
the experimental range of stability of the rock salt crys- 
tal structure. The groundstate energies of three different 
magnetic structures: (i) ferromagnetic (F), (ii) type II 
antiferromagnetic and (iii) antiferromagnetic with prop- 
agation vector (0,0,^), were calculated self-consistently 
and mapped on Heisenberg Hamiltonian 



H — 2_^ Jij ^i ' ^j 



(1) 



with nearest-neighbor Ji and next-nearest-neighbor J2 
interaction, which is known experimentally to provide a 
good description of magnetic behavior of the materials 
in question. The classical energies corresponding to spin 
configurations (i) to (iii) are — (12Ji+6J2)5'^, 6J2S^ and 
(4Ji — 6J2)6'^ respectively. The ferromagnetic energies 
were calculated with the lattices of both antiferromag- 
netic structures and the corresponding energy differences 
were taken. The assumption of J1 — J2 model was checked 
for EuS at ambient pressure by spin-spiral calculation 
similar to those of Ref. [l^, which yielded dispersion con- 
sistent with Ji — J2 model. The calculated exchange 
constants as function of lattice parameter for values of U 
from 6 eV to 9 eV are shown in Fig. Q and Fig. [3 

From the parameters of the Heisenberg Hamiltonian 
ordering temperatures were calculated using result of the 
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FIG. 2: The same as in Fig. [Hfor EuSe and EuTe. 



Tyablikov decoupling method |20|, l2l| 



where J(q) stands for the lattice Fourier transform of 
the effective exchange parameter. The alignment in the 
groundstate is determined by the sign of Ji+ J2 (positive- 
ferromagnetic, negative-antiferromagnetic). For EuTe 
we have calculated the Neel temperature as well using 
generalized equation ((2j 22]. The ordering temperatures 
as functions of lattice constant are shown in Fig. [31 and 
Fig. ^ and compared to the experimental data of Gon- 
charenko and Mirebeau [23, 12^ (see also |3: 14)- The 
estimated numerical accuracy is about 10 K. The trend 
of weakening ferromagnetism in favor of antiferromag- 
netism when going from oxide to telluride is well repro- 
duced and is common to all values of U. The effective 
exchange is rather sensitive to the value of U, yet the 
literature value of 6 eV to 7 eV 25] gives the best agree- 
ment throughout the series. Both the significant increase 
of the Curie temperature with pressure and almost con- 
stant behavior of the Neel temperature of EuTe at lower 
pressures observed experimentally are well captured by 
the calculations 

Based on the fact that in the type II antiferromagnetic 
structure the first neighbor exchange is frustrated and 
thus the mean-field Neel temperature is proportional to 
J2 , Goncharenko and Mirebeu concluded that J2 is pres- 
sure independent while Ji exhibits a non-linear increase 
with the applied pressure. In their scenario the transi- 
tion from low-pressure antiferromagnetic groundstate to 
high-pressure ferromagnetic groundstate is solely due to 
the increase of Ji. Our calculations provide a different 
picture. With the only exception being EuO at high pres- 
sure, we find more or less linear dependence of Ji on the 
lattice parameter. On the other hand J2 exhibits quite 
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FIG. 3: The magnetic ordering temperatures of EuO and 
EuS calculated for different values of U (see caption of Fig. 
QJ. The experimental values for EuO ,2^ and EuS |2J are 
marked with open diamonds. 



non-linear behavior, which in the case of EuTe translates 
to being almost constant at low pressures and increas- 
ing rapidly at higher pressures, which significantly con- 
tributes to stabilization of the ferromagnetic state. 

The success of the LDA-hU functional in describing the 
trend across the chalcogenide series as well as capturing 
the pressure dependence of the ordering temperature in- 
dicates that the relevant coupling mechanisms are well 
accounted for. While we can not identify directly the 
leading coupling mechanisms we can make several obser- 
vations connected to the sensitivity of the results to the 
value of U. The value of U affects inter-site coupling in 
two distinct ways: (i) through the position of the occu- 
pied / bands within the semiconducting gap, which is de- 
cisive for hybridization with the valence and conduction 
bands, (ii) through the splitting between the occupied 
and unoccupied / bands. In order to stress the impor- 
tance of including hybridization effects in the / bands we 
show the spin-majority bandstructure of EuO in Fig. [SJ 
The / band exhibits appreciable dispersion with a band- 
width of about 1 eV which is mostly due to mixing with 
O 2p valence bands. A similar picture is obtained for 
the other members of the series with the / bands being 
localized deeper in the gap when going toward telluride. 

The intra-atomic f — s or f — d exchange is determined 
solely by the spin density. The coupling mechanisms in- 
volving only the intra-atomic exchange are therefore in- 
sensitive to the value of U. The Bloembergen-Rowland 
(BR) coupling [23, inter-band analog of the Ruderman- 
Kittel-Kasuya-Yosida (RKKY) coupling well known in 
metals, is the leading mechanism of this type. A dis- 
cussion of BR coupling was given earlier by Lee and 
Liu p , starting from diagonalized band Hamiltonian and 
adding interband exchange term, and by Kasuya |2i] and 
Sawatzky et al. [^, starting from atomic orbitals with 
f — d exchange and conduction- valence (Eu d - X p) hy- 
bridization. The sizable dependence of our results on 
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FIG. 4: The magnetic ordering temperature of EuSe and 
EuTe calculated for different values of U (see caption of Fig. 
QJ. The negative values for EuTe correspond to Neel temper- 
ature and antiferromagnetic groundstate. The experimental 
values marked with open diamonds are taken from Ref. [2^. 




FIG. 5: Majority spin bandstructure of EuO obtained with 
U of 7 eV at ambient pressure. The valence bands have dom- 
inant O 2p character while the conduction bands are mostly 
Eu d and s states. 



the value of U indicates that additional mechanisms are 
involved. 

We have performed a simple model calculation to eval- 
uate the role of unoccupied / bands which participate, for 
example, in the superexchange mechanism. By adding an 
orbital dependent auxiliary potential, which acts only on 
the minority-spin / orbitals, we control the 'site energy' 
of the unoccupied (minority-spin) / states without effect- 
ing the occupied /'s. Obviously such a term does not en- 
ter the groundstate energy directly, but only through the 
mixing of the minority-spin / bands with the occupied 
bands. Excitations from occupied to unoccupied bands 
correspond to hopping of an / electron from one Eu atom 
to another, resulting in superexchange interaction. By 
looking at the dependence of the effective exchange pa- 
rameters on the splitting of the occupied and unoccupied 
/ bands controlled by the auxiliary potential we can as- 
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FIG. 6: Difference between the groundstate energies of type 
II antiferromagnetic and ferromagnetic structures as a func- 
tion of the sphtting Uf between the occupied and unoccupied 
/ bands obtained with auxihary potential on top of LDA+U 
with U of 7 eV. Uf of 11 eV corresponds to zero auxihary 
potential. 



sess the relative importance of the superexchange and 
virtual processes involving change in the 4/ occupation 
in general. In Fig. El we show the energy difference be- 
tween the ferromagnetic and type II antiferromagnetic 
groundstates as a function of energy separation between 
the occupied and unoccupied / bands. Apparently the 
auxiliary potential has a sizable effect consistent with 
l/Uf dependency of the superexchange interaction. 

Now we summarize. Our calculations show that ac- 
counting for intra-atomic repulsion using the LDA+U 
method provides a reliable description of effective ex- 
change coupling in ferromagnetic insulators with lo- 
calized moments. The trend favoring ferromagnetism 
for lighter chalcogenides as well as the under-pressure 
antiferro-to-ferromagnetic transition in EuSe and EuTe 
are well captured. The pressure dependences of the mag- 
netic ordering temperatures, which correspond well to 
the experimental observations, are connected to distinct 
under-pressure behavior of the exchange parameters Jl 
and J2. The exchange coupling is strongly effected even 
by the unoccupied / bands establishing the important 
role of virtual hopping processes between the / shells on 
neighboring atoms. 
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